Two Phase Collective Modes in Josephson Vortex Lattice in Intrinsic Josephson 

Junction E^S^CaC^Og.^ 
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Josephson plasma excitations in the high T c superconductor Bi2Sr2CaCu2 0s+« have been investi- 
gated in a wide microwave frequency region (9.8 - 75 GHz), in particular, in magnetic field applied 
parallel to the ab plane of the single crystal. In sharp contrast to the case for magnetic fields par- 
allel to the c axis or tilted from the ab plane, it was found that there are two kinds of resonance 
modes, which are split in energy and possess two distinctly different magnetic field dependences. 
One always lies higher in energy than the other and has a shallow minimum at about 0.8 kOe, then 
increases linearly with magnetic field. On the other hand, another mode begins to appear only in a 
magnetic field (from a few kOe and higher) and has a weakly decreasing tendency with increasing 
magnetic field. By comparing with a recent theoretical model the higher energy mode can naturally 
be attributed to the Josephson plasma resonance mode propagating along the primitive reciprocal 
lattice vector of the Josephson vortex lattice, whereas the lower frequency mode is assigned to the 
novel phase collective mode of the Josephson vortex lattice, which has never been observed before. 



I. INTRODUCTION 

The phase collective excitation in Josephson junctions 
known as Josephson plasma has been studied since the 
1960's, triggered by a theoretical prediction of Ander- 
soni However, there are only a few experimental re- 
ports on Josephson plasma resonance in single junctions 
made by conventional superconductors^ because there 
is a strong damping effect of the quasiparticles due to 
small superconducting energy gap in the single junc- 
tions. This situation was improved drastically in the 
case of high T c superconductors because they are com- 
prised of well-established weakly coupled Josephson junc- 
tions along the crystallographic c axis (intrinsic Joseph- 
son junction), whose plasma frequency along the c axis is 
much lower than the superconducting gap. In particular, 
in the case of Bi 2 Sr 2 CaCu20 8+5 (Bi2212), 4 ' 5 the super- 
conducting CuC>2 layers and the insulating or semicon- 
ducting Bi202 layers are regularly stacked in an atomic 
level in the unit cell forming a weakly coupled Josephson 
junctions. Based on this, Tachiki and co-workers pre- 
dicted the existence of the Josephson plasma excitation 
in high T c superconductors*!^ Subsequently, two theoret- 
ical models have been developed : one by Koyama and 
MachidafSai the other by Pedersen and Sakai4S*ii Quite 
recently, a unified theory has been reported*^ 

The Josephson plasma resonance in Bi2212 is observed 
as a sharp and strong microwave absorption in a finite 
magnetic field parallel to the c axis^ Especially, the 
longitudinal Josephson plasma mode gives an extremely 
sharp resonance because of its small dispersion and the 
high quasiparticle damping rate of Bi2212 14-15 . The 
clear observation of the longitudinal plasma resonance 
has been a strong advantage of Bi2212 for investigation 
of the Josephson plasma phenomena in supercon- 
ductors compared with many other materials such 
as Las-^Sr^CuOzti 16 Bi 2 (Sr,La)2CuO 6+ , 5 (Bi2201) f iI 
YBaaCuaOr-^YBCO)^ and BEDT-TTF salts^ in 
which the Josephson plasma resonance (edge) is less 



clear. 

According to the theory of the Josephson plasma reso- 
nance, the resonance frequency uj p (H,T) can be written 
asLSa 



w 2 p (H,T) = ujI (T) (cos m+1 (H,T)) : 



(1) 



where to p {T) = c/y / eA c (T) is the Josephson plasma 
frequency in the absence of a magnetic field. e 
and A C (T) stand for the dielectric constant and the 
temperature-dependent c axis penetration depth, respec- 
tively. ipi^i+i(H,T) is the gauge-invariant phase differ- 
ence between the (Z)th and (/ + l)th layers and (• ■ ■ ) 
denotes the spatial and time averages. 

In the case of a magnetic field parallel to the c axis 
(perpendicular to the superconducting CuC>2 layers), 
it is well known that pancake vortices are generated 
above H c \, the lower critical field. For the Josephson 
plasma resonance, this situation is taken into account 
by considering the uniform reduction of the Josephson 
current along the c-axis, resulting in the decrease of 
(cos(pij+i(H, T)), — i.e., the reduction of U) P (H,T) in 
Eq. QJ. Because of this fact, the Josephson plasma 
resonance can be used as a sensitive method to evaluate 
(cosy M+1 (g,T)), 21 i 22 i 23 . 24 

This treatment cannot be applied to the case where the 
external magnetic field is parallel to the ab plane. In this 
case, the Josephson vortices are introduced in-between 
the CuC>2 double layers and contribute the Josephson 
plasma resonance as an essential ingredient of the phe- 
nomenon. In contrast to the case for the perpendicular 
field, the plasma frequency is hardly suppressed because 
the fluctuations of pancake vortices are too small to dom- 
inate the interlayer coupling. This is clearly observed in 
the angular dependence of the Josephson plasma reso- 
nance near the ab plane. A sharp dip in the Josephson 
plasma resonance field is clearly observed within an an- 
gle of ± 2 degrees from the ab planed From this sharp 
change of the Josephson plasma resonance behavior, it is 



2 



in turn expected that the dynamical motion of Joseph- 
son vortices such as Josephson vortex lattice collective 
modes has been thought to account for the Josephson 
plasma resonance frequency. This opens entirely new 
possibilities for the Josephson plasma resonance phenom- 
ena in the case of parallel magnetic field configurations 
and deserves more detailed consideration both experi- 
mentally and theoretically. Since the Josephson vortex 
modulates the interlayer Josephson tunnelling current in 
the length scale of Xj = js along the ab plane, where 
7 is the anisotropy parameter and s is the interlayer 
distance, the collective plasma oscillation strongly cou- 
ples with the collective Josephson vortex motion such as 
Josephson vortex lattice modes. 

This picture is initiated by the single junction model, 
which was studied in the 1960s by Lebwohl and Stephen^ 
and by Fetter and Stephen^! theoretically. They solved a 
one- dimensional nonlinear wave equation and obtained 
two modes as excitation spectra: one of them lies above 
the plasma gap huj p in zero field and approaches asymp- 
totically the linear relation with the propagation vector 
k (plasma mode) . Another mode obtained by Fetter and 
Stephen in a single junction is a gapless vortex sliding 
mode. The plasma mode with an array of Josephson vor- 
tices was first of all observed indirectly as a so-called Eck 
resonance in the current -voltage (I—V) characteristics of 
the Pb/PbO/Pb single junction, in which the resonance 
was detected at a voltage proportional to the applied 
parallel field H\\~ Fiske has also shown an anomalous 
step-like behavior in the I — V characteristics in Al, Sn, 
Pb, and Nb junctions^ and this phenomenon was asso- 
ciated with the resonant electromagnetic modes of the 
junction. 

More recently, Josephson plasma resonance experi- 
ments have been performed in the external magnetic field 
configuration being closely parallel to the ab plane in 
single-crystal Bi2212iS^2i2i They observed anomalous 
behaviors such as the angular dependence of the reso- 
nance field associated with the sudden jump of it, which 
have partly been explained by the theory of Bulaevskii 
et alM> It seems that the resonance mode indicating this 
anomalous behavior can be described qualitatively by the 
single junction modei^l. However, the model in which 
conventional (metallic) single junctions are connected in 
series cannot be applied in intrinsic Josephson junctions 
especially in treating Josephson plasma excitations be- 
cause the charge conservation inside a layer does not hold 
due to thinner CuC>2 double layers (3 A) than the charge 
screening length (fi - 10 A). Therefore, we consider that 
the previous treatment of the Josephson plasma reso- 
nance in a parallel magnetic field is insufficient and ex- 
clude the rich physics inherent to the intrinsic Josephson 
junctions. 

In this paper, we report observation of two-phase col- 
lective modes through systematic measurements of mi- 
crowave absorption as functions of magnetic fields paral- 
lel to the ab plane, temperature, and frequency in Bi2212 
single crystals, and then the origins of the two modes 



are argued. The organization is as follows. Experimen- 
tal results by using the microwave resonance technique 
described in Sec. [H]are presented in Sec. IIIII In this sec- 
tion, we first present the experimental data at 25.5 GHz, 
where the phenomena are the most interesting and the 
richest among all frequencies (Sec. 1111 Af) . Next, we show 
the experimental results at different frequencies and de- 
scribe the frequency dependence in order to obtain a field 
dependence of the resonance modes (Sec. IIII B|) . Finally, 
universal features deduced from the doping dependence 
are described (Sec. IIII Cll . The origins of the observed 
two modes are discussed in Sec. II VI with the help of re- 
cent theoretical modelsi^^i^ 



II. EXPERIMENTAL SETUP 

In order to perform Josephson plasma experiments in 
the Josephson vortex state, it is necessary to have a 
wide range of microwave frequencies beyond the zero-field 
plasma frequency, which is expected above 100 GHz for 
optimally doped Bi2212. Since in our experimental facil- 
ity the microwave frequencies were limited up to 90 GHz, 
it was necessary to bring down the zero-field Josephson 
plasma resonance in the range of 50 GHz. This was 
achieved by shifting the doping level to the underdoped 
region by annealing the pristine overdoped Bi2212 sam- 
ples under reduced atmospheres. 

We have measured microwave absorption in three 
Bi2212 single crystals grown by the modified traveling 
solvent floating zone (TSFZ) method. Two underdoped 
crystals (Ul, U2) and one optimallydoped crystal (OP) 
were used in the measurements for the comparison. The 
superconducting transition temperatures T c 's of these 
crystals were determined by low-field magnetization mea- 
surements with a superconducting quantum interference 
device (SQUID) magnetometer as 70.2, 76.8, and 90.5 K 
with superconducting transition widths of 1.0, 2.5, and 
0.8 K for Ul, U2, and OP, respectively. The typical size 
of the crystals is 0.8 x 0.8 x 0.02 mm 3 . 

Most of the microwave absorption data except for those 
in Sec. IIII Gl have been obtained in Ul. Measurements 
were made in a frequency range between 9.8 and 75 GHz 
by both reflection- and transmission-type bridge bal- 
ance techniques using rectangular cavity resonators with 
TE102 mode. In order to excite the Josephson plasma, 
the samples were placed inside the cavity in such a way 
that the oscillating electric field of the microwaves was 
exerted parallel to the c axis over the ab plane of the crys- 
tal as shown in Fig. 1 of Ref. 36 . Frequency-stabilized mi- 
crowaves were generated by the synthesized signal gener- 
ator (Hewlett Packard 83650B) or Gunn oscillators, and 
the magnetic field was applied by a split-pair supercon- 
ducting magnet. The angle between field direction and 
Cu02 plane of the sample was adjusted by rotating the 
cavity resonator with respect to the field direction by a 
precision rotator within an accuracy of 0.001 degree. 

The resonance data were obtained either by sweeping 
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FIG. 1: The resonance curves obtained at 25.5 GHz by sweep- 
ing temperature (a) and by sweeping field (b). In (a), temper- 
ature was swept from far above T c to the lowest temperature 
after changing magnetic field at each measurement, and no 
hysteresis was found. In (b), the solid and broken curves in- 
dicate data obtained in increasing and decreasing fields. Hys- 
teresis found in field sweep measurement becomes more signif- 
icant at lower temperatures and in lower fields. Thick arrows 
in (b) denote shoulders attributed to the weak resonance of 
HTM above the turnover. The sharp peak shown in (b) at 10 
kOe is due to DPPH as a field marker. 
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FIG. 2: The resonance curves observed at 9.8 (a), 18.9 (b), 
and 22.3 (c) GHz. The temperature is swept down from far 
beyond T c after setting magnetic fields. At 9.8 GHz (a), the 
separation of LTM and HTM cannot be identified although 
the line-shape changes asymmetric to symmetric at a low field 
region below 0.5 kOe. At higher frequencies than 18.9 GHz, 
a sharp HTM and relatively broad LTM are split well and 
easily identified. 



III. EXPERIMENTAL RESULTS 



A. Simultaneous observation of two resonance 
modes 



the magnetic field at various fixed temperatures or by 
sweeping the temperature at various fixed fields. The res- 
onance occurs when lo p (H,T) matches with the incident 
microwave frequency u> by varying either the magnetic 
field or temperature. It is worth noting that the field 
sweep (FS) measurements would give a non-equilibrium 
Josephson vortex state because the strong hysteretic be- 
havior of the resonance was observed, especially at low 
temperatures while the magnetic field was swept up and 
down. In order to avoid this unnecessary complication, 
we employed measurements mostly by the temperature 
sweep (TS) when low-temperature data were needed. 
The temperature was either stabilized within 20 mK dur- 
ing the field sweep or swept at a rate of 2 K/min. It is 
also worth noting that since it is technically difficult to 
sweep microwave frequency as far as the cavity resonator 
method is employed, we are compelled to prepare a num- 
ber of cavities (about 15 cavities) with different resonance 
frequencies. 



Figure^a) shows microwave absorption curves at 25.5 
GHz in various constant parallel magnetic fields H\\ ob- 
tained by TS measurements. Two clear resonance lines 
with different characters are found at higher and lower 
temperatures. As is seen clearly, there is a well-separated 
temperature gap in-between. This feature strongly de- 
pends on the microwave frequencies: the higher in fre- 
quency the larger separation in temperature is observed. 
As will be shown later, the lower temperature mode 
(LTM) quickly shifts towards the low temperature side 
and disappears with increasing microwave frequency w. 
On the other hand, with lowering frequency below 30 
GHz the LTM quickly shifts up and merges with the 
higher temperature mode (HTM), which is explained in 
Fig. |21 In general, the HTM in this frequency region ap- 
pears only in a relatively low field, approximately below 
3 kOe, as long as u> is lower than 55 GHz, which is close 
to the zero-field and zero-temperature plasma frequency 
u>p(0) described later. When ui goes beyond 33 GHz, only 
one resonance is observed as shown in Fig. [3| In the fre- 
quency region between 10 and 30 GHz both resonances 
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FIG. 3: The resonance curves observed at 34.5 (a), 52.3 
(b), and 74.8 (c) GHz by sweeping temperature from above 
T c . In this high frequency region, only HTM was observed. 
The resonance drastically shifts to lower temperatures and 
broadens as u is increased. At 74.8 GHz, HTM was obtained 
only above 3 kOe without zero field resonance. 



come into a play, resulting in a complicated phenomenon 
as a function of temperature as well as magnetic field. 

Focusing on the HTM, it behaves in a very unusual 
manner as seen in Fig. ^a), for example. In zero mag- 
netic field a sharp resonance with a symmetric lineshape 
with respect to temperature was clearly observed. Since 
this resonance corresponds to the zero-field mode argued 
previously, the temperature dependence of the zero-filed 
plasma frequency lo p (T) can be described by the two- fluid 
model assuming the conventional Drude models 

In a finite parallel magnetic field below approximately 
1 kOe, HTM begins to shift slightly to lower tempera- 
tures. With further increase of the magnetic field, how- 
ever, it turns to shift backward to higher temperatures 
at around 0.8 kOe as obviously seen in Fig. ^a) [this 
behavior is more clearly displayed in Fig. 0] (b), where 
the data obtained at many frequencies measured are pre- 
sented for the LTM]. The strong absorption intensity of 
this HTM in zero field quickly decreases as the field is 
applied. The lincwidth in terms of temperature is wider 
in finite fields than in the zero field. 

It has been well-established that uj p is always 
suppressed by application of perpendicular magnetic 
fields?2222 because the Josephson current due to coher- 
ence effect between layers is a decreasing function of mag- 
netic fieldiSi This requires that in higher magnetic fields 
the resonance temperature must be lower in order to gain 
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FIG. 4: (Color online) (a): The resonance field-temperature 
plot for the LTM resonance obtained at 9.8, 18.9, 22.3, and 
25.5 GHz. (b): The resonance field-temperature plot for the 
HTM resonance obtained at 9.8, 18.9, 22.3, 25.5, 34.5, 39.5, 
44.2, 52.3, 61.7, 65.9, 74.3, and 74.8 GHz from right to left. 
The same symbols represent data obtained at the same fre- 
quencies. Inset: The lo v vs. T plot at zero field. Asterisks are 
experimental data and the solid line is a curve fitted to the 
two fluid model. 



the fraction of Josephson current which is reduced by the 
perpendicular magnetic fields. It seems that this is not 
the case in the Josephson plasma resonance in parallel 
fields to the layers and this phenomenon certainly re- 
quires a new explanation. 

Such an unusual turnover behavior just mentioned 
above cannot be explained either without making an un- 
usual assumption that the phase coherence between lay- 
ers would increase further with further increase of the 
parallel magnetic field approximately above 0.8 kOe or 
without introducing entirely new mechanisms for the 
Josephson plasma resonance in the parallel magnetic 
field. It is worth mentioning that this increasing ten- 
dency of the temperature dependence of H ICS toward 
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FIG. 5: The resonance curves at 9.8 (a), 18.9 (b), and 22.3 
(c) GHz by sweeping magnetic field to H = after stabi- 
lizing temperatures at the maximum fields. The resonance 
peaks pointed by inverted triangles correspond to the HTM 
resonance. 



higher temperatures seems never to go beyond T c in zero 
field, although it goes beyond the plasma resonance tem- 
perature To in zero field at the frequency u> as seen in 
Fig. Etb). It fades away quickly just below T c . This sur- 
prising result may allow us to speculate that the HTM 
is more stable in temperatures with the formation of the 
Josephson lattice state than without Josephson lattices. 
This unusual behavior has never been observed in the 
Josephson plasma resonance in perpendicular magnetic 
fields studied earlier and is certainly a new phenomenon 
only observable in parallel magnetic fields. 

In parallel magnetic fields, consecutive occupation of 
the Josephson vortices in-between the Cu02 layers is ex- 
pected to be realized along the c axis as a function of 
increasing magnetic field, preserving an isosceles triangu- 
lar symmetry according to the theoretical calculations. 38 
This formation of Josephson vortex lattices would cause 
the oscillatory coherent Josephson current, which accord- 
ingly possesses the wave vector k associated with the pe- 
riodicity of the Josephson vortex, and would provide the 
resonance at even higher temperatures in a finite field 
than that in zero field according to the periodic Joseph- 
son vortex arrangement. 

At the intermediate field region between 1.3 and 2.2 
kOe both the HTM and LTM can be observed at the 
same magnetic field as seen in Fig. [TJa). For instance, 
the data at 1.69 kOe indicate the existence of both the 
HTM and LTM with 25.5 GHz at (1.69 kOe, 66.2 K) 
and (1.69 kOe, 43.0 K), respectively. This feature is ob- 



served only within a certain frequency region between 18 
and 30 GHz, and only one resonance line can be found 
at the other frequencies. This observation of two reso- 
nance modes indicates a peculiarity in the the Josephson 
plasma resonance in H || ab unlike the one for H || c, 
where only one resonance line as far as the longitudinal 
mode is concerned, was observed to monotonically de- 
crease with both increasing temperature and magnetic 
field. 

The resonance curves obtained by FS measurements 
are displayed in Fig. H(b). As the temperature is just 
below T c (= 70.2 K) a sharp resonance begins to appear 
(To ~ 68 K) from the zero field limit. The resonance 
rapidly grows with decreasing temperature. As a few 
degrees Kelvin below To, two resonance line can be dis- 
tinguished: one at a lower field is stronger in intensity 
than the other one which lies at a higher fields and at 
the foot of the resonance of the stronger line (pointed by 
arrows). This two-resonance feature can only be found 
in 62 < T < 66 K, below which the resonance suddenly 
changes its character with broadening and with highly 
the hysteretic nature as seen in Fig. ^b). The temper- 
ature range is wider at higher frequencies. For example, 
it extends to 14 K at 52.3 GHz, above which frequency 
the two resonance feature can no longer be observed. 

As seen in Fig. QJb) the absorption curves show a 
considerable hysteretic behavior for the field being swept 
up and down, especially below 45 K. This hysteretic fea- 
ture is not appreciable at high temperatures above 60 K 
whereas it is most prominent in the intermediate temper- 
ature region between 10 and 45 K. It is noticed that there 
are several small resonancelike peaks at the lower-field 
side in addition to the main large resonance. However, 
this anomaly as well as even in the main peaks depends 
strongly on the field sweeping direction and the sweep- 
ing rate. Furthermore, it must be mentioned that some 
faint peaks are observed in-between the HTM and LTM 
in the FS mode [at 52 K in Fig. ^ (b), for example], 
where there is no resonance in the TS mode because of 
the gap between the HTM and LTM as seen in Fig. Q] 
(a). The peaks are, perhaps, non-resonant absorptions 
due to the non-equilibrium effect of the Josephson vor- 
tex state because the peak fields are almost equal to the 
fields where the LTM is observed at the highest temper- 
atures. Therefore, we do not consider these faint peaks 
in the FS measurement hereafter. 



B. Frequency Dependence 

1. Low-frequency region (uj/2tt < 30 GHz) 

The typical resonance curves at 9.8, 18.9, and 22.3 GHz 
obtained by TS measurements are shown in Figs. |2Ia), 
|5{b), and^c), respectively. At the lowest frequency of 
9.8 GHz, only one resonance line is observed. The reso- 
nance line shape changes from asymmetric to symmetric 
with increasing field below 0.5 kOe, indicating a precur- 
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sor effect of splitting into the HTM and LTM already at 
this frequency region. In sharp contrast to this, the reso- 
nance splitting is clearly observed at 18.9 GHz as shown 
in Fig. |2Ib), where the LTM resonance and the HTM res- 
onance appear below and above 64 K, respectively. The 
HTM resonance becomes weak in intensity and is quickly 
smeared with increasing field, and finally disappears at 
1.0 kOe, while the LTM resonance appears from the low- 
temperature side above about 0.7 kOe, then shifts to the 
higher-temperature side, and then turns back to lower 
temperatures with increasing field with little change in 
the lincshape up to 3 kOe. The larger separation between 
the LTM and HTM resonances at higher frequencies indi- 
cates that the LTM has a different frequency dependence 
from the HTM. The integrated H — T diagrams for most 
frequencies for both the LTM and HTM are shown in Fig. 
21 by plotting the resonance temperatures as a function 
of field. 

The resonance curves obtained in FS measurements for 
various temperatures at the same frequencies shown in 
Fig. |2Jare shown in Fig. |SJ All resonance peaks displayed 
here belong the LTM except for the ones indicated by 
the inverted triangles. The resonance peaks of the LTM 
in increasing field lie at a lower field than the ones in 
decreasing field, and the hysteresis observed in increasing 
and decreasing magnetic fields becomes more significant 
at lower temperatures and in lower fields. In comparison 
with the TS data, the resonance peaks obtained in the 
case of decreasing field coincide with the peaks in the TS 
measurements rather well, in which more homogeneous 
Josephson vortex state is expected to occur. 

It is very unlikely that a complete lock-in vortex state 
is realized without any pancake vortices in the actual 
experimental condition, because the calculated angle to 
be allowed to form the lock-in state by taking into ac- 
count of the sample geometry (the interlayer distance of 
12 Aand the length of the sample of 0.8 mm) would be 
of the order of 10~ 4 deg or less, which is much smaller 
than the experimental angular resolution. Furthermore, 
from the x-ray parallel-beam double-crystal method, it 
has been known that our sample used for the present ex- 
periment has a width of the rocking curve of the order 
of 0.025°. This means that the crystallographic imper- 
fectness of our sample, especially the parallelism of the 
Cu02 layers, gives an actual limitation to the measure- 
ments. Therefore, we think that the Josephson vortex 
state is always accompanied by some amount of the re- 
maining pancake vortices. In such a situation the vortex 
crossing lattice state would be formed as first suggested 
by KoshelevfSi where the interaction between Josephson 
vortices and pancake vortices may be attractive. Since 
the pancake vortices are pinned strongly at low temper- 
atures as seen in the large hysteresis in magnetization, 
the Josephson vortex (JV) system would be compelled 
to be pinned by the pancake vortices through an attrac- 
tive interaction. At the present stage of knowledge we 
do not know to what extent this pinning effect is impor- 
tant in the Josephson plasma phenomena except for the 



hysteretic effect while the magnetic field is swept. 

However, it is important to note that the complete JV 
state can be established in a smaller sample with the di- 
mensions less than 10/im x 10/im in the ab plane. This 
is evidenced by the angular dependence measurement of 
the JV flow resistance which shows a sudden rise at an 
angle within one degree or so from the ab plane, depend- 
ing strongly on the sample size and the intensity of the 
applied magnetic field ~ 



2. High frequency region 

Figures |3f a) , |3f b) , and |3£ c) depict resonance absorp- 
tion curves at 34.5, 52.3, and 74.8 GHz, respectively, by 
TS measurements at various fixed magnetic fields. In 
this higher-frequency region, only the HTM is observed. 
As seen from Fig. [31 it is indispensable to note that 
the HTM at 34.5 and 52.3 GHz has zero field resonance, 
while it does not have zero field resonance at 74.8 GHz. 
This is because the HTM in the low-field region quickly 
shifts to the low temperature side and disappears above 
approximately 57 GHz. This can be seen more clearly in 
Fig. Efb), where the HTM is shown in the H -T plot. 
This zero-field resonance occurs only at a particular tem- 
perature To at the corresponding to. The uj dependence 
of To has been reported previously^ 41 - and is explained 
by using a simple two-fluid model with a temperature- 
independent scattering rate^ By extrapolating this tem- 
perature dependence of the plasma frequency uj p {T) to 
T = as shown in the inset of Fig. Efb), the inherent 
(zero temperature) plasma frequency lu p (0)/2tt is esti- 
mated to be 56.8 GHz, which yields A c = 217/im with 
e = 15. This value of A c is quite reasonable for under- 
doped BSCCO with T c of 70.2 K and is used to quantify 
the field dependence of the HTM in the later sections. 

In the frequency range below but near u/ p (0), the dra- 
matic turnover behavior of the resonance temperature is 
observed as seen in Fig. Efb) • The Josephson plasma res- 
onance at low fields below 1 kOe shifts to the lower tem- 
perature side and then turns back to higher temperatures 
above 1 kOe as the magnetic field is increased. A typical 
example is shown in the H—T plot in the left panel of Fig. 
0in the case of uj/2tt = 44.2 GHz, where T = 55.7 K and 
the shift in temperature amounts to about 3.5 K below To 
at about 1 kOe, then the resonance peak sharply turns 
over toward the high-temperature side, and reaches to 
59 K, which is even more than 3 K higher than Tq. The 
intensity of the resonance becomes weaker and weaker 
with increasing magnetic field, and the resonance finally 
fades away. As the frequency is increased, the whole 
curves including To moves towards low temperature as 
shown in the H — T plot in Fig. 0fb). However, it is 
interesting to note that the disappearing temperature of 
the resonance after the turnover lies in the temperature 
region between 55 and 65 K even above the zero- field 
resonance temperature Tq. This surprising behavior is 
very different from our understanding of the Josephson 
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field-temperature plot shown in Fig. |3]to the ui-H diagram 
at various frequencies at 55 K as an example. 



plasma resonance for H || c, where the plasma frequency 
should decrease with application of a magnetic field and 
with increasing temperature as a result of the suppres- 
sion of the interlayer coherence (cos<pij+\{H,T)) due to 
the fluctuation of pancake vortices as formulated in Eq. 

At frequencies above u> p (0), there is no longer the res- 
onance line at H = but only one resonance in a fi- 
nite field in the HTM. As shown in Fig. 01c), the reso- 
nance line can be observed only above about 3 kOe, and 
the resonance occurs at higher temperatures in higher 
magnetic field until reaching the limiting temperature of 
about 60 K around 5 kOe. This means that the excitation 
energy of the HTM becomes higher at higher tempera- 
tures. Since the superconducting phase coherence should 
become weaker and weaker at higher temperatures, this 
increasing excitation energy of the Josephson plasma in 
the HTM needs a new expression if it is enhanced fur- 
ther by overcoming the reducing temperature effect. As 
shown in Fig. 01b), three more data sets obtained at 
different frequencies are included. 



3. Overall features 

In order to construct the u> p — H diagram, which is a 
direct indication of the H\\ effect on the excitation modes, 
the following procedure is applied: Since the most of ex- 
periments have been done in the TS mode at a fixed 
magnetic field at a certain u> and the H — T diagram was 
constructed as shown in Fig. one should convert this to 
the tUp — H plot by taking the resonance points from cor- 
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responding to at various temperatures as shown in Fig. 
|SJ The resonance frequency is plotted as a function of 
magnetic field for nine different temperatures in Fig. 
It is noted that the plots in Fig. 0were extracted from 
resonance field versus temperature plots at all measured 
frequencies (partly not shown in Fig. 0} . 

It is clear that there are two different type of resonance 
modes with very different field dependent characteristic 
features: one corresponds to the HTM, which starts with 
u>p(0) at zero field; then, the frequency grows monoton- 
ically and almost linearly after passing a shallow mini- 
mum in the uj p vs H diagram as the field is increased. 
On the contrary to this, the LTM shows a weak initial 
increasing behavior with increasing magnetic fields and 
then turns to the gradual decreasing behavior after pass- 
ing a broad maximum. This LTM seems to be disappear 
at the zero-field limit. 

As mentioned above, the HTM deviates from linear 
behavior and tends to have a minimum at 0.5 - 1.0 kOc, 
and then merges with the zero field Josephson plasma 
mode. Since this HTM is directly connected to the zero- 
field Josephson plasma mode, the HTM is confirmed to 
be the Josephson plasma mode in the Josephson vortex 
state in a parallel external magnetic field. It is noted 
that the initial decrease of w p (H\\) is stronger at higher 
temperatures. On the other hand, the LTM exhibits a 
broad peak at around 1.5 kOe and then tends to decrease 
again with decreasing field. It is important to note that 
this LTM disappears at the limit of zero magnetic field. 
This experimental fact implies strongly that the existence 
of JV's is crucial to the LTM. 

In Fig. [SI u> p (Hn) normalized by the zero field plasma 
frequency at a given temperature derived from the inset 
of Fig. 01 b) is shown for five typical temperatures in Fig. 
As seen in Fig. |S]it is interesting to point out that 
both the HTM and LTM are scaled well to two single 
curves by using us p {T). 



C. Universality of Two Excitation Modes 

The microwave measurements have been performed in 
a similar manner in other crystals with higher doping 
levels (U2 and OP samples). The essential features de- 
scribed above were unchanged: the characteristic two 
modes HTM and LTM exist with qualitatively similar 
temperature as well as magnetic field dependencies. A 
typical example of the experimental data are shown in 
Figs. EJa) and E{b) for an underdoped U2 sample and 
for an optimally doped OP sample, respectively. Com- 
paring these data with the data shown in Figs. [2] and [31 
one can easily find that the separation between the HTM 
and LTM for U2 is smaller than that of Ul at a similar 
frequency. For example, the data in Fig. EJa) (at 34.4 
GHz) can be compared with the data in Fig. EJa) (at 
34.5 GHz) and are rather similar to the those obtained 
at 22.3 GHz shown in Fig. Efc). Assuming that these 
two data sets obtained at 22.3 and 34.5 GHz provide the 
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FIG. 9: The resonance data obtained (a) at 34.6 GHz in 
U2 (T c =76.8 K) and (b) at 44.1 GHz in OP (T c =91.0 K). 
Two resonances and the temperature gap are found in both 
samples. Separation of two modes at 4 kOe in U2 (a) is similar 
as the data at 22.3 GHz and 1.6 kOe in Ul (Fig. H (c)). 
The inset of (a) shows temperature dependence of the zero- 
field plasma frequency in U2. Dotted lines trace the HTM 
resonance and is guide for the eyes, which shows turnover 
behavior as well as HF mode in Ul. 



same ratio of ui/u> p (0) to each sample, u> p (0)/2tt for U2 
can be derived as 88.1 GHz. This estimated value of 
u) p (0) is consistent and agrees rather well with the value 
extrapolated by the temperature dependence of the zero- 
field Josephson plasma frequency uj p (T) as shown in the 
inset of Fig. H^a), although the number of experimental 
points are only a few. 

In the case of the OP sample, typical experimental 
data are shown in Fig. Etb) at 44.1 GHz, where no clear 
splitting of HTM and LTM is observed. Applying the 
similar scaling in this case (for 18.9 GHz in Ul), a value 
of ui p (0)/2ir w 130 GHz is derived. This value is close 
to the one obtained for an optimally doped Bi2212 single 
crystal4i Thus we suggest that the resonance splitting 
can be noticeably observed at 10 > w p (0)/3 and the mag- 
nitude of the T-gap is larger for u>/u> p (0) being closer to 
unity. In contrast, u> is much smaller than tv p (0), and 
the LTM merges to HTM as observed at 9.8 GHz in Ul, 
resulting in a single branch of the H — T plot. 

In previous studies, Matsuda et alm& and Tsui et alM- 
have reported that the Josephson plasma resonance field 
has a sharp decrease near 9 = (H || ab). In particu- 
lar, Matsuda et al. found the disappearance of resonance 
in an underdoped crystal (T c = 87 K) through measure- 
ments of the angular dependence of the FS data and ar- 
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FIG. 10: A schematic picture of propagation of the tilted 
plasma wave. The interlayer phase difference along the ab 
plane and the c axis are spatially modified as shown by dotted 
and broken waves, respectively. The propagation vector of the 
plasma wave comprised of the transverse component q x and 
the longitudinal component q z is depicted as a thick arrow. 



gued the existence of a collective mode of the Joseph- 
son vortex lattice which lies above the experimental mi- 
crowave frequency (u/2ir — 45 GHz) at zero field and 
increases with parallel magnetic fields on the basis of a 
theory by Bulaevskii et al£&. Judging from our results, 
this vanishing of the resonance observed by Matsuda et 
al. can be interpreted differently because their sample 
is underdoped crystal and the measurement was done at 
36 K. Their experimental conditions seem to be in the 
gap region between the HTM and LTM. Furthermore, 
according to their results the disappearance was not ob- 
served in an optimally doped crystal (T c = 89.5 K). This 
also seems to be consistent with the systematics which we 
found here although they attributed it to the misalign- 
ment of the magnetic field. From an experimental point 
of view, no matter what the theoretical interpretation 
may be, the observed dramatic change of the Josephson 
plasma resonance near 9 — is strongly suggestive for the 
excitation of the new phase collective modes in intrinsic 
Josephson junctions. 



IV. ORIGIN OF EXCITATION MODES 

A. High Temperature Mode 

As discussed by Fetter and Stephen in single 
junctions, 27 propagating Josephson plasma waves are 
strongly modified by the JV's, resulting in the Josephson 
plasma mode having a linear field dependence according 
to the linear increase of the reciprocal lattice vector of the 
JV array, kn = 2irsH\\/$o, where <£>o is the flux quan- 
tum. This indicates that a periodic arrangement of JVs is 



crucial to have the propagating Josephson plasma mode. 
We first assume that JV's penetrate all block layers and 
form isosceles triangular lattice with the lattice constants 
along the c axis and the ab plane being 2s and $o/si/|| 
in high-i?|| region, respectively, as shown in Fig. 1101 al- 
though for layered superconductors the argument has not 
been completely settled yet. A theoretical calculation by 
Ichioka 42 predicts that such a dense JV lattice can be 
realized above H* = \/3^q/8^s 2 , which gives 1.8 kOe 
for 7 = 1085 obtained from A c = 217/xm and X a t — 2000 
A. This value is rather in good agreement with the mag- 
netic field, where the HTM begins to tend to the linear 
relation and the LTM has a broad maximum. Hence, it 
is reasonable to think that the HTM above H* may be 
the Josephson plasma resonance mode in the dense JV 
lattice. 

Considering a Josephson plasma wave traveling in the 
isosceles triangular JV lattice, the k vector of the Joseph- 
son plasma mode corresponds to the primitive reciprocal 
lattice vector of the JV lattice q, which consists of the 
constant c-axis component q z = tt/s and the H» propor- 
tional a6-plane component q x = 2ttsH^/<5> , as depicted 
in Fig. Here, y and z axes are taken parallel to the 
JV and c axis, respectively, and x is taken perpendicu- 
lar to the y and z axes. Thus, Josephson plasma waves 
can be excited at all possible wave vectors as of the mix- 
ture of the longitudinal plasma with k z = q z and the 
transverse plasma with k x = q x . The plasma frequency 
is determined by the dispersion relation along q lying 
between the longitudinal and transverse plasma modes 
and increases with magnetic field because of the linear 
increase in q x . It should be noted that the dispersion is 
very close to the longitudinal dispersion for q x <C q z . 

Bulaevskii et al. have proposed a formulation of the 
plasma resonance in a layered superconductor on the ba- 
sis of the single- junction modeljF 2 - and Koshelev and 
Machida have derived excitation spectra by consider- 
ing the longitudinal coupling effect In the high field 
limit, where all block layers are occupied by JV's form- 
ing an isosceles triangular lattice, the peak frequency of 
the dissipation spectrum, corresponding to the resonance 
peak, as a function of parallel magnetic field is formulated 
as 



Wp(-ff||) tt^s 2 

w P (0) ~ 



(2) 



Using 7 = 1085, this equation provides an excellent 
agreement with the experimental data above 3 kOe at 
all temperatures, as shown in Fig. [5J It is also consistent 
that the linear field dependence of the HTM is violated 
below H*. We therefore conclude that the linear increase 
in Wp(H\\) of the HTM is due to an increase in q a b of the 
JV lattice, while q c is fixed on tt/s because of the intrin- 
sic pinning. This result suggests that the excited plasma 
frequency can be easily controlled by adjusting the mag- 
netic fields parallel to the layers even above lo p . 
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B. Low Temperature Mode 

The LTM lies at uj < ui p /2 and slightly decreases with 
Uii as displayed in Fig. [H] The LTM in crystals with 
higher doping tends to appear at higher frequencies al- 
though the whole u> p — H diagram has not been estab- 
lished yet. In under doped YBCO, a similar excitation 
mode was observed around uj p /2 with slight negative Hu 
dependence and is referred to the a mode^. The a mode 
cannot be observed after zero-field cooling, suggesting 
the considerable influence of pinning of the JV's on the 
collective oscillation mode. This may also be the case 
in Bi2212 as shown in Fig. ^b) and as well as in Fig. 
|SJ where a large hysteretic effect is observed. Assuming 
that our LTM corresponds to their a mode, the low-lying 
collective mode should strongly depend on u> p , which is 
a material parameter and given by the strength of the 
Josephson coupling between CuC>2 bilayers. 

It is well known that the lowest collective mode of 
a Josephson vortex system is the vortex sliding mode 
with k x = 0i2i The sliding mode becomes gapless in an 
ideal system without pinning and couples with a dc ho- 
mogenous electric field perpendicular to the supercon- 
ducting layers. Since the pinning of Josephson vortices 
cannot completely be removed in actual samples, the low- 
est mode should have a nonzero frequency. In a junction 
with a periodic pinning for instance, the Josephson crit- 
ical current depends periodically on x£l 

Since the hysteresis obtained by the FS measurements 
indicates that the JV pinning strongly affects to the 
LTM, we once discussed the pinned vortex sliding sce- 
nario as a possible origin of the LTM. 44 However, our 
finding that LTM tends to become higher in more highly 
doped samples with higher lu p cannot be explained by 
this scenario, because it is not natural that the collective 
mode hardened by the pinning effect can be scaled by w p . 
Therefore, the vortex sliding model is not suited for and 
not likely for the origin of LTM in the framework of con- 
ventional knowledge, and needs for a new explanation. 

In order to reveal the origin of the LTM, it is re- 
quired to introduce an intrinsic stacking effect of intrinsic 
Josephson junctions because previous theoretical models 
giving the vortex sliding mode are based on the single 
junction model, where the correlation of the JV lattice 
along the z axis is neglected. The vortex sliding mode in 
layered superconductors which we have in mind means 
coherent vortex motion parallel and perpendicular to the 
layers with k x = k z = 0. Here, we may consider a collec- 
tive mode with k x — but k z ^ at a finite frequency. 
Considering an antiphase oscillation of JV arrays between 
adjacent junctions with k z — n/s and k x — 0, for in- 
stance, the oscillation is a sharing motion between adja- 
cent JV arrays as a result of dynamic phase oscillations 
not only between adjacent layers but also next-adjacent 
layers, so that the frequency of the oscillation can be 
considered to be scaled by io p because the interaction be- 



tween JV's in adjacent junctions is to be governed by the 
Josephson coupling as the first approximation. 

Very recently, Koyama reported theoretical calcula- 
tions on the Josephson plasma resonance in the JV 
latticed After numerical calculations on analytically 
derived equations, a collective mode which lies about 
u p (0)/2 in the low-field limit and slightly decreases with 
increasing field was obtained. He argues that the collec- 
tive mode is considered as an anti-phase oscillation mode 
of JV's which originates from the strong charge coupling 
between junctions, and it is also responsible for the lon- 
gitudinal Josephson plasma mode. This picture seems 
to be most appropriate as a candidate for the origin of 
the LTM although no direct experimental proof has been 
obtained. 

For deeper understandings of the LTM, further theo- 
retical and experimental investigations are needed in the 
low-frequency region. We think that an investigation of 
the JV flow branch in the current-voltage characteristics 
along the c axis is a complementary experiment to eluci- 
date the JV behavior in the layered superconductors in 
the low-frequency region. 



V. CONCLUSIONS 

We have clearly identified two microwave excitation 
modes with a temperature-dependent gap in parallel 
magnetic fields in single crystal Bi 2 Sr 2 CaCu208+5. It 
is shown that the strong coupling effect between the 
Josephson plasma and JV lattice is responsible for these 
two modes. The high frequency mode is attributed to the 
Josephson plasma mode propagating along the primitive 
reciprocal lattice vector of the JV lattice. Making use of 
this interplay, we may be able to explore electromagnetic 
waves in the frequency range between the superconduct- 
ing gap and ui p . The low frequency mode is attributed to 
the new collective mode which has never been observed 
before and seems to be unique for layered superconduc- 
tors. It is considered that the Josephson vortex collective 
oscillation mode with finite k z contributes to the phase 
oscillations, although a quantitative full understanding 
has not been obtained. 
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